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Abstract 

The spin of a free electron is stable but its position is not. Recent quan- 
tum information research by G. Svetlichny, J. Tolar, and G. Chadzitaskos 
have shown that the Feynman position path integral can be mathemati- 
cally defined as a product of incompatible states; that is, as a product of 
mutually unbiased bases (MUBs). Since the more common use of MUBs 
is in finite dimensional Ifilbert spaces, this raises the question "what hap- 
pens when spin path integrals are computed over products of MUBs?" 
Such an assumption makes spin no longer stable. We show that the usual 
spin- 1/2 is obtained in the long-time limit in three orthogonal solutions 
that we associate with the three elementary particle generations. We give 
applications to the masses of the elementary leptons. 

The first section discusses mutually unbiased bases and position path in- 
tegrals, and the difference in behavior between position and spin. Section 2 
introduces spin path integrals over MUBs. Section 3 derives some convenient 
arithmetic results for products of spin-1/2 projection operators. Section 4 cal- 
culates the long-time MUB spin path integrals. Section 5 shows that the long 
time propagators, after summing over orientation, converge to the usual spin- 
1/2. Section 6 applies the results to the lepton masses. Finally, section 7 
discusses the results. 

1 Introduction 

Let A = {\aj)} and B — {\bk)} be two different bases for a finite dimensional 
Hilbert space. They are "mutually unbiased" if all inner products {aj\bk), have 
the same magnitude. If the dimension of the Hilbert space is N, then the inner 
products have magnitude y^l/N : 



\{a,\h)\ = Vl/N. 



(1) 



The simplest example of mutually unbiased bases are spin- 1/2 in two perpen- 
dicular directions. Bases for spin-1/2 in the x, y and z direction are: 




These three bases are mutually unbiased; the magnitudes of the transition am- 
plitudes are all y/T/N = ■\/l/2. This is a complete set; there are only three 
perpendicular directions in three dimensions. 

Zee's textbook introduction to quantum field theory [ij introduces the path 
integral formulation with the double slit experiment. "A particle emitted from 
a source S at time t — passes through one or the other of two holes, Ai and 
A2 drilled in a screen and is detected at time t = T by a detector located at O. 
The amplitude for detection is given by a fundamental postulate of quantum 
mechanics, the superposition principle, as the sum of the amplitude for the 
particle to propagate from the source S through the hole Ai and then onward 
to the point O and the amplitude for the particle to propagate from the source 
S through the hole A2 and then onward to the point O." 

Increasing the number of holes increases the number of paths. Thus with 
three holes the total amplitude is the sum of three single path amplitudes: 

A{S -^0) = ^^iA{S Aj O). (3) 

Adding two more screens B and C, between A and O, see Fig. ([iJ, requires 
summing the amplitudes over 3'^ = 27 paths: 

A{S ^ O) = I]|^iI]LiSLiyl(S' -> Aj Bk -^Ci^ O). (4) 

The path integral formalism follows by considering increases in the number of 
intermediate screens and the number of holes drilled in them. If the screens 
have enough holes, and are sufficiently closely spaced, one requires that the 
calculation should give the transition amplitude for a free particle moving in 
empty space between S and O. 

For a particle moving in a potential V{q), the final result for the amplitude 

is 

{qT\e-'"^\qo) = J Dqit)e' -I'o -""^^^l (5) 

In the above Dq(t) represents all possible paths q which run from the source 
S = g(0) to the detector O = q{T). The contribution from each path consists 
of a complex phase. 

In 2008, G. Svetlichny 0] noticed that while MUBs and path integrals appear 
to be very different things, they both involve amplitudes where the information 
is contained only in the phases. He showed that in the limit of short time 
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Figure 1: Paths from 5 to O through three intermediate screens A, B, and C, 
with three holes each. 

intervals, the path integral approaches the transition amplitudes between two 
MUBs. Thus longer paths consist of products of MUB transition amplitudes. 
In 2009, J. Tolar and G. Chadzitaskos [3] obtained the free particle propagator 
as a sum over products of MUB transition amplitudes by a limiting procedure 
in finite dimensional Hilbert spaces. 

To see quantum behavior in the motion of an electron we must measure its 
position with an accuracy smaller than its de Broglie wavelength: 



For an electron with energy 1 keV, this distance is A = 0.4 A = 4 x 10~^^m. Of 
the particle detectors we have available, the most accurate is emulsion which 
can measure particle positions to an accuracy of around 5 x lO^^m. [Hj This is 4 
orders of magnitude larger than the de Broglie wavelength of a 1 keV electron. 
Heavier and higher energy particles have even smaller de Broglie wavelengths. 
Consequently, elementary particle tracks appear classical. @ Instead, the best 
evidence we have for the bizarre behavior of quantum particles over short times 
and accurate positions is obtained from diffraction experiments such as the 
single slit and double slit experiments. 

In the very short time limit, a product of MUB transition amplitudes ap- 
proaches a single MUB transition amplitude. In this case all possible transitions 
are equally probable. This behavior corresponds to the familiar result of single 
slit experiments: with a sufficiently narrow slit, the particles receive random 
velocities and their tracks spread out. 

The free propagator for a spin-1/2 particle does not change spin; that is, the 
transition amplitudes between spin- up | -I- 2), and spin-down | — z) are zero. In 
terms of path integrals over spin space, the paths do not cross, spin-up stays 
spin- up. See Fig. If we measure the spin of a beam of spin- up electrons the 
result is always spin-up. This paper considers the possibility that if we could 
measure spin over a sufficiently short time interval, we would find the same 
behavior as position: the measurement of spin would modify the spin. Thus 
the traditional Stern-Gerlach measurement of spin-1/2 is classical in the sense 



A = 



12.2 



A, 



(6) 
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that continuous particle tracks are classical. We assume that the underlying 
quantum behavior consists of transitions between mutually unbiased bases. For 
an interesting argument that the Stern- Gerlach experiments can be interpreted 
entirely from a classical understanding of electricity and magnetism, see 0, 
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Figure 2: A free particle beginning with spin-up stays that way; forbidden 
transitions shown with sparse dots. 

The effort here is similar (and uses similar mathematics) to that of Foster and 
Jacobson; see and references therein. They attempt to obtain the massive 
Dirac propagator from the chiral left and right handed states. In doing this, 
they consider paths that move along one of four directions in spacetime. The 
steps are oriented towards the corners of a tetrahedron in 3 dimensions and move 
one unit forwards in time. Since they are replacing the usual two-element basis 
set for spin-1/2 with a set of size four, one would expect them to obtain two 
copies of the Dirac propagator. They avoid the unwanted doubling by making 
computations that begin and end with spinors. Thus their propagators amount 
to 2 X 2 matrices and describe only the desired two states. These ideas originated 
with Feynman's "checkerboard" model of the electron in 1 -|- 1 dimensions, [loj . 

The present paper's calculation begins with a tripled basis set for spin-1/2 
and since we're looking for three generations of massless chiral fcrmions, we do 
not avoid the tripling. Instead of spinors, our calculations concentrate on what 
happens to the propagators at long time. In addition, we work in the quantum 
information approximation; we do not consider position. 

2 Spin Projection Operators 

For the Hilbert space with dimension M=2, Svetlichny's prescription for the 
analog to the Feynman path integral is (see equation (9) of Q): 

{qrlqa) = ^xi^X2--'^XM{<lT\xi){xi\x2)-{xN-l\xN){xN\qo) (7) 

where we require that xj and Xj+i be mutually unbiased for all j. This sum 
over products reduces since for each j, we have that the set {xj} are a complete 
set: Y^x:j\xj){xj \ = 1. 

Equation ^ is unsatisfactory in that it implies a fixed sequence of choices 
of basis state. For example, the first sum Sa;^ , might be over spin in the ±x 
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direction, the next sum over spin in the ±z, etc. While it may be mathematicaUy 
true, it is not physically satisfying. 

If we eliminated the arbitrary choice of basis by allowing Xj to run over any 
of the six states in the three MUBs given in Eq. © we would eliminate the 
arbitrary choice of bases, but the calculation would still give a trivial result 
since Yia:-\xj){xj \ — 3. Instead we will require that Xj be taken from one of the 
three positive spin states: 

l+x), 1+2/), \+z). (8) 

Thus instead of summing over the two states taken from one of the three (ar- 
bitrarily chosen) basis states, we will allow our sums to run over three basis 
states, one chosen from each basis. 

We will be summing quantities such as the following: 

{+x\ + ^{+z\+x){+x\+y){+y\+z). (9) 

Starting on the right, this corresponds to a particle whose spin path goes through 
the sequence | + z) to | + j/) , to | + x) , to \+ z), and finally to | + x) . To make 
these calculations, we will rewrite the above as: 

{+x\ \ + ^{+z\ I + x) {+x\ I + V) (+y I \+z). (10) 

The operator |+a;)(+a;| is the projection operator for spin in the +x direction 
and similarly for +y and +z. We will label these spin projection operators X, 
F, and Z. In terms of the Pauli spin matrices, they can be written as: 

X = |+f)(+f| = (l + c7,)/2, 

Y = \+y){+y\ = (l + ^y)/2, (11) 
z = \ + ^{+A = (i + ff.)/2. 

In this notation, Eq. (O becomes {+x\XZXYZ\ + z). 

More generally, for any vector u, define cr„ = Ux(Tx + UyGy + u^cr^. Then 
(1 + <Ty)/2 is the projection operator for spin-1/2 in the +u direction, and 
(1 — CTti)/2 projects spin-1/2 in the — u direction. We will use a bar for the 
opposite spin state, C7 = (1 — cr„)/2. 

When one converts a basis set {| + u), | — u)} into projection operators, the 
features which characterize a basis set are transformed into algebraic relations 
among the projection operators. The normality of the kets become idempotency: 

{±u\ ± u) = 1 ^ [(1 ± a„)/2] [(1 ± (7„)/2] = (1 ± a„)/2, (12) 

and orthogonality becomes annihilation: 

{±u\ T m) - ^ [(1 ± a„)/2] [(1 T ^n)/2] = 0. (13) 

Finally, the requirement that the number of basis elements equals the dimension- 
ality of the vector space becomes the requirement that the sum of the projection 
operators is unity: 

(l-fa„)/2 + (l-r7„)/2 = l. (14) 
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In addition the projection operators have trace 1: 

tr((l ± n^)/2) = tr(l/2) ± tr(a„/2) = 1. (15) 

We wih find the same relations among the long-time MUB spin path integrals. 
Our set of annihilating projection operators that sum to unity will have three 
basis elements; we will associate these with the three generations of elementary 
spin- 1/2 fermions. 



3 Spin Projection Arithmetic 

There are three products of three spin projection operators that begin and end 
with Z: 

Z X Z = y/lp2?Z, 

ZY Z = VT722Z, (16) 

z z z = \ftj¥z. 

These are real multiples of Z . The powers of two in the above equations count 
the number of transitions. With ZXZ and ZYZ there are two transitions, and 
with ZZZ = Z there are none. More general products can be complex: 



ZXY Z= Z, (17) 

where we define ^/S = exp(±z7r/4). A path that makes N transitions through 
different projection operators will have a magnitude of V2~^ . 

The complex phase in Eq. (ITTl) is a geometric phase, also called Berry 



ll| or Pancharatnam [12| phase. These phases can be picked up when a quan- 
tum particle goes through a series of states and returns to its initial state. It 
does not depend on the arbitrary complex phases of spinors and consequently is 
an observable. Our use of quantum phase will arise from products of projection 



operators of the sort described in R. Bhandari's paper |13|, but for spin-1/2 
rather than photon polarization. 

In computing sums of MUB Feynman paths, we will use X, Y, and Z as 
part of a basis for a complex vector space. These three account for paths that 
begin and end with the same projection operator. Paths whose final projection 
operator are different from their initial projection operator must be handled 
differently. For these mixed paths there are six ways to choose different initial 
and final projection operators. 

There are three possible products of three projection operators that begin 
with Z (on the right) and end with X: 



X X Z = 072" XZ, 

XYZ = ^/+i/¥XZ, (18) 
X Z Z ^ XZ. 

In the above calculation, the three products are all complex multiples of the 
same matrix: XZ . Twice this matrix is idempotent and has trace 1: 

i2XZr = 2XZ= do). (19) 
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This is a property of all MUB paths that begin and end with different projection 
operators; they will be complex multiples of an idempotent with unit trace. 

Products that begin and end with orthogonal projection operators are the 
familiar raising and lowering operators. For example ZXZ and ZY Z are both 
raising operators for spin-1/2. They differ by a complex phase: 

ZXZ = e^'I'^ZYZ. (20) 

The path ZXZ ZY Z picks up the same complex phase, 7r/2: 

ZXZYZ = e"/2 ^1/24 Z. (21) 

Thus the choice of phase in raising and lowering operators are an example of 
geometric phase. 

The pure density matrices are usually defined as the matrices that can be 
produced from normalized spinors. An alternative definition is that they are the 
Hermitian projection operators with trace 1. To convert a pure density matrix to 
a spinor choose a nonzero column (or the complex conjugate of a nonzero row), 
treat the elements as a vector, and normalize. The matrix used in Eq. ([T9| . and 
the other five obtained by replacing X and Z with different projection operators, 
possesses two of the three properties that define a pure density matrix; they are 
idempotent and have trace 1 but they are not Hermitian; their rows give different 
spinors than their columns. 

Hermiticity is associated with the property of time reversal invariance. We 
have: 

{2XZ)^ = 2ZX. (22) 

More generally, Hermitian conjugation reverses the order of a product of any 
number of projection operators. This corresponds to the reversal of the time 
ordering of a path. 

For convenience, we will drop off the factor of 2 and use 

{X, y, Z, XY, YX, XZ, ZX, YZ, ZY} (23) 

as the basis for MUB spin path integrals. Any product of MUB projection 
operators can be written as a complex multiple of one of these nine. In using 
these, we have five more than are needed to give a basis for 2 x 2 complex 
matrices, but what we are looking for is a basis for paths and path integral 
calculations. 

In Feynman path integrals, one specifies initial and final states of the par- 
ticles and sums over diagrams that connect those states. For the MUB case 
we will do the same. For this reason, our set of nine products of projection 
operators can be used as if they were a basis set for a complex vector space 
of dimension nine. A vector in that space is a collection of nine amplitudes. 
Adding two such vectors together is equivalent to adding nine Feynman path 
integrals to nine other Feynman path integrals to get nine sums. In this way, 
given a path, (or a collection of MUB paths that all share the same initial and 
final projection operators), we associate a complex number. 
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In addition to summing two paths, we also need to connect one path to 
another. This is equivalent to the multiplication of two products of projection 
operators, for example, {XY){YZ) = (XYYZ). Suppose we have two paths (or 
sums of paths) A and B which are associated with complex numbers a and b. If 
it so happens that path A ends with the same projection operator that path B 
begins with, then the concatenation of these paths, call it BA, will be associated 
with the complex number ba = ah multiplied by a non commutative correction. 
The correction depends only on the initial and final projection operators of the 
paths. This suggests that we should organize our calculations so that transitions 
between projection operators are kept inside products of the path integral basis 
Eq. I|23|) : we can then use complex multiplication to model the concatenation 
of paths. 

Accordingly, we will only concatenate two path integral basis elements if 
the final projection operator of the first element matches the initial projection 
operator of the second element. For example, this requirement allows {XY){YZ) 
or {Z){ZX) but not {XY){Z) or {Z){X). Since the path integral basis has nine 
elements, there are a total of 81 possible products of them but the restriction 
reduces the number of products we will consider to 27. Below, we will show 
that this is sufficient to account for any spin path, see Eq. ([55)) . 

Of the 27 remaining products of path integral basis elements, 15 are already 
correct for standard complex multiplication. The five that begin with X are: 

{X){X) - X, 
{X){XY) = XY, 

{X){XZ) = XZ, (24) 
iXY){Y) = XY, 
{XZ){Z) = XZ. 

The other ten such products are obtained by cyclic permutation of X, Y, and 
Z. We will call these the "diagonal products" for reasons which will be clear 
below. The remaining 12, "off diagonal products" are more complicated. The 
four that begin with X are: 



{XY){YX) = VW^, 

{XZ){ZX) = y ^X , 

{XY){YZ) = y /+i/¥ XZ, 

{XZ){ZY) = y/^t/¥ XY. 



(25) 



The other eight are obtained by cyclic permutation. 

Other than the complex factors of the off diagonal products, the above 27 
products are compatible with matrix multiplication. Let {a'^,a'y, ...a'y^) be a 
9-vector of complex numbers associated with a collection of MUB spin paths. 
We use the prime to indicate that these represent complex multiples of non 
commutative matrices. For the moment we will suppose that there is some 
single path associated, perhaps of type XZ, so that only one of them (i.e. a'^^) 
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is nonzero. Assemble them into a 3 x 3 matrix a' as follows: 



(26) 



In this matrix, the columns represent the three initial states (+j7|, {+z\, 

while the three rows represent the final states | + a;) , \ +y), and \+ z). 

Let b' be a similar matrix for another path whose final projection operator 
matches the initial projection operator of a', that is, Z. For instance, b' could 
represent a path of type ZY so that only b'^y is nonzero. Concatenating the 
paths gives a path of type XY. Corresponding to this, the matrix product (ab)' 
will have only one nonzero entry, {aby^y — a'^zKy- This is not quite correct; 
according to Eq. ([25]) we should have {ab)'^y — a'^^b'^y. 

To fix the matrix product, we need to scale the off diagonal elements (i.e. 



' „l n' „l n' 



^Tlll ^^1^1 



xy^ yz^ ^zx^ ^yx^ zy 



a' , a'^^ and similarly for b') in such a way that the off diag- 



onal products are corrected without changing the diagonal products. Rewriting 
Eq. (|25l) in terms of what it says about products of the elements of a' and 6', 
what we want are new matrices, a, b, and (ab) such that: 



(27) 



71/2^ (a6); = a'^yb'y,, 

Vl/2^{aby, = aLKx, 

V+W{aby^z = <yb'yz^ 

^J~^iJ^ {aby^y = a'^zKy^ 

becomes {ab)x — axybyx^ etc. A suitable transformation is a' — > a by: 

a-x = a^, 

axy = % a'^y, (28) 
* / 

O-xz — fjg 0-XZ1 

where 

rjg = VT72 e+'^/i^e^^f'^/s, for g = 1, 2, 3. (29) 

Cyclic permutations give the transformations on the other six path basis ele- 
ments. The complex phase 2i7r/3 appears repeatedly; we will abbreviate it as 
w: 

■w = exp(2i7r/3), (30) 

so ?7g = \/l/2 exp(i7r/12) . We will use the integer parameter g to represent 
the generation quantum number. 

The matrix entries of b' represent the nine cases: spin changing from -fa; to 
+x, +x to +y, etc. Suppose this propagator is followed by another propagator 
a' . In computing the propagator (a6)' we must sum over all possible paths. With 
matrices a and b, this simply amounts to matrix multiplication. For example, 
the three terms on the right side of 

{ab)xz = axbxz + axybyz + axzbz- (31) 
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correspond to the three paths that go through the new node with spin 

and +z, respectively. Thus we have transformed the problem of concatenating 

the propagators of MUB spin into 3x3 matrix multiplication. 



4 Long-Time MUB Spin Propagator 

Let's begin by computing the spin path integral (or propagator) from +x to +z 

with two internal paths, call it 6*2 (+a;, +z). This is a sum over paths with four 
projection operators. The initial projection operator is X, the final projection 
operator Z while the two inner projection operators can be any of X, Y, or Z. 
Taking all possible cases for the inner projection operators and calculating with 
the Pauli spin matrices, or using the rules given above, we find: 

G2{+x,+z) = ZXXX + ZXYX + ZXZX 
+ ZYXX + ZYYX + ZYZX 
+ ZZXX + ZZYX + ZZZX, 



(3 - 3j/4) 



1 l^ (^2) 




= (6 - 3i/2)ZX. 

This amplitude is too large to conserve probability. 

The problem is that the transition amplitudes need to be adjusted for the 
fact that we have added two new possible paths at each vertex. To preserve 
probability at each vertex, we need to make the transition probabilities smaller. 
There are three transitions so we will multiply by a factor k^. 

We will compute k later in this section, for now, let's see how to rewrite a 
path as a matrix multiplication. We are concerned with paths that look like 
{Z){P){Q){X) where (P) and (Q) can be any of the three projection operators. 
First, using idempotency, we duplicate P and Q and factor into pairs: 

ZPQX = ZiPP)iQQ)X, 

= {ZP){PQ){QX). ^""""^ 

Modify the pairs to replace (XX) = (X), {YY) = [Y), and [ZZ) = (Z). Now 
we have the path as a product of our nine path basis elements in such a way that 
each pair of adjacent basis elements match their adjacent projection operators. 

The state ZP has final state \ +z), and initial state anything, so we represent 
it as K times the matrix with Is in the bottom {z) row: 




[zp)' = K 1^ I . (34) 
Convert the matrix {zp)' to {zp) to obtain: 

(zp) = k\ 1 (35) 
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Similarly, the middle term {PQ) represents any initial state and any final state 
so it will be a matrix (pq)' with all elements equal to k. This is the general 
propagator for a single MUB step. Transforming it to (pq) we have: 



(pq) = k\ 7]* 1 T]g \ ■ (36) 




Finally, (QX) will be k times a matrix whose left (x) column only is nonzero, 
with all entries 1. It converts to: 




(qx) \ v; , (37) 

V % ; 

and the sum over paths is represented by the matrix product: 

G2(+.f,+i) =kM I I r;* 1 % I I r;* I . (38) 




^5 

Upon multiplying, and factoring out the ZX matrix to the right, we obtain: 

/ \ 

G2(+f,+f) = St,^[i%+V*gV*g+V*gr]^g)/Vg]\ 0. (39) 

V % J 

Replacing jjg with ■\/l/2exp(i7r/12)z«3, we find that 

^K^^iVg + rfgV; + V*gVg)/Vg = «'(6 " 3*/2). (40) 

The result does not depend on g and is the same as the sum over paths, Eq. ([5^ . 

The above result generalizes. Longer path integrals introduce extra factors 
of the (pq) array of Eq. ([55)) and extra factors of k. Define: 

1 Vg V*g 

- « I 1 % I , (41) 



Vg Vg 



1 



so longer paths will involve powers of Gg. Note that Gg is 1-circulant, that is 
each row is identical to the row above, but rotated 1 position to the right. 

The discrete Fourier transform diagonalizes 1-circulant matrices. This allows 
us to compute G^ by taking the discrete Fourier transform, taking powers of 
the diagonal entries, and then reverse transforming. Thus we can solve for G^ 
in closed form. 

Define the Fourier transform matrix F as 



(42) 
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so that the discrete Fourier transform of a vector v is Fv. Then the discrete 
Fourier transform of a matrix M is 

M = FMF*. (43) 

Given a 1-circulant matrix with top row {A, B, C), the discrete Fourier transform 
converts it to a diagonal matrix: 

A + w*B + wC 

A + wB + w*C I . (44) 

A+B+C 

The elements down the diagonal are \/3 times the discrete Fourier transform of 
the vector (C, B, A)*. The transform of Gg is a diagonal matrix Gg of this form. 
The jth element on the diagonal of Gg is: 

[Gg]jj = + \/2cos(257r/3 - 2jir/3 + 7r/12)]. (45) 

The three matrices Gg have the same diagonal elements but in an order that 
depends on g. 

The largest entry on the diagonal of Gg is [Ggjgg] it will dominate (Gg)^ as 
— > oo. In order for the limit to exist, this diagonal entry must be 1. Therefore 
we have that 

1 = [Gg]gg = + V2cos(7r/12)]. (46) 

The other two diagonal entries in Gg are [Gg]g+i^g+i = k(1 + •\/2cos(— 77r/12)) 
= 2 - \/3 and [Gg]g+2,g+2 = k[1 + N/2cos(-157r/12)] = 0. We have: 

/I 

G^=F* {2-V3)^ ) F, (47) 
\ 

and similarly for and G^, but with the elements on the diagonal rotated. 
Finally, to obtain the non commutative amplitudes, one performs the reverse 
transformation. 

In the limit as iV goes to infinity, G^ becomes a diagonal matrix with the 
^th entry equal to one and all other entries zero. Taking the inverse discrete 
Fourier transform we find: 



1 w-s 
y+9 1 ur- 
yj-g y^+a 1 



G~ = F*(Gg)°°F = - I w+3 1 w-3 |. (48) 



The magnitude of all the entries are 1/3 so the transition probabilities all equal 
(1/3)^ = 1/9. On converting the above to non commutative form, we find: 



V2 e-'^'/^^w- 



Gg°° = ^ I e+*"/i2y^+9 y/2 e-'^'/^^wS | . (49) 
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The above can be translated from the matrix using the path basis to obtain: 



G'g°° = [{X + Y + Z) + e-'"l^'^w-3.JYj2(XY + YZ + ZX) ^^^^ 



+e+"/i2y^+9yT72(rx + ZY + xzy^i-i. 

We wiU associate these three propagators with the three generations of elemen- 
tary fermions. 



5 The Particle Generations 

The long-time propagators of Eq. (|50l) are idempotent: 

(g;-)2 = G-;-, (51) 

they annihilate each other: 

G;°°Gr = 0, iig^K (52) 

and have unit traces. The off diagonal elements of G^ are each multiplied by 
a factor of w"*"^ or . Since the sum of powers of w is zero, i.e, 

w + w"^ +w'^ =w^^ +w^'^ +w^^ =Q, (53) 

the three matrices Gg°° sum to the unit matrix: 

G'{° +Gi" +G^ =X + Y + Z = 1. (54) 

These are the same relationship Eq. (IT^ to Eq. (ITSI) , exhibited by the projection 
operators for a complete set of basis states. Therefore, ignoring the short-time 
behavior, the long-time propagators define a 3-dimensional Hilbert space. The 
three projection operators X, Y , and Z define another 3-dimensional Hilbert 
space. This is just three times the number of states in the usual spin-1/2. We 
associate the tripling with the three generations. 

The factors and that distinguish the three generations are cubed 
roots of unity. Since the long-time propagators are otherwise identical, this 
suggests that the generations should be simple when characterized as cube roots 
of unity. That is, w"*"^ are the three solutions to the equation: 

= 1. (55) 

Better, the off diagonal elements of Eq. are 

Ug = exp(2i57r/3 -I- i7r/12)/V2. (56) 

These are the roots of the equation z'^ — (l-|-i)/4. We expect that the differences 
between the generations of elementary particles should be simple when expressed 
as functions of Ua- 
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What we have done amounts to complexifying generation; we extend gener- 
ation from a discrete variable that takes on the three values 1, 2, and 3, to a 
complex variable that can take on any complex value. This extension of gen- 
eration amounts to consideration of particle states that are not the orthogonal 
long-time propagators. This is similar to the work of T. Regge who extended 
orbital angular momentum to a complex variable. 

Our long time propagators contain extra information beyond what is present 
in a spin-1/2 state. Instead of a spin axis, we are provided instead with a set of 
three orthogonal unit vectors x, y, and z. From these we obtain the spin axis 
as the normalized sum; 

u^{x + y + z)/V3. (57) 
The projection operator for spin in this direction u = (1, 1, 1)/a/3 is 

c7u = il + {a^+ay+a,)/V3)/2. (58) 

Our problem is that other complete sets of MUBs will give the same spin axis 
and spin projection operator. 

Given two spin- 1/2 states \+u) and \+v), the transition probability between 
them is given by the square of the inner product: 

\{+u\ +v)\^ =tr(| +u){+ u\\ + v){+v\\ + u){+u\). (59) 

The right hand side of the above is written in terms of matrix multiplication of 
pure density matrices. Our long time propagators Eq. (j48l) are written in matrix 
form and their concatenation is given by matrix multiplication so we can use 
this form to verify that the long time propagators lead to the usual spin- 1/2 
relations. 

In computing Eq. ()59|) . we need to make sure that the rows and columns of 
our matrices are compatible in terms of their orientations. In general, | -I- u) 
and I -|- v) will not share the same complete MUB sets so we cannot repeat our 
trick, used in Eq. (1331) . of arranging for the matrix multiplication to be between 
identical bases. Instead, we will have to sum over orientations. 

MUBs must remain perpendicular; consequently there are only three orien- 
tations possible - those corresponding to the even permutations of the three 
bases: 

{xyz^yzx, zxy}. (60) 

The nine complex numbers in the long time propagators are left unchanged by 
these permutations. Consequently, when we sum over these permutations, the 
result will be alterations of the path basis Eq. ([23l) . For example, the three 
orientations will take the path basis element X to {X, Y, Z}, so the result of the 
sum will heX^X + Y + Z, which can be rewritten as {3 + ax + (Ty + az)/'2. 
Similarly, Y and Z will be taken to X + Y + Z. The complete transformation 
on the path basis is: 

X,Y,Z ^ X + Y + Z = {3 + iax+ay+a,))/2, 
XY,YZ,ZX XY + YZ + ZX = {'S + (2 + i){a^ + dy + az))/A, 

YX,ZY,XZ YX + ZY + XZ = {3 + {2-i){a^+(7y+(7,))/A. 

(61) 
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The right hand sides of the above are hnear combinations of 1 and {(7x+(7y + az) 
which is close to what we need, Eq. ([58|. 

Summing over orientations will reduce the 3x3 matrix form of a long time 
propagator to a 1 x 1 form, that is, it will be reduced to an object that can be 
represented among the usual complex 2x2 matrices. The result Eq. (|5T|) shows 
that this object will be some complex linear combination of 1 and au- 

Before we can compute the spin-1/2 interaction between two long time prop- 
agators we need to renormalize the long time propagators for spin-1/2 type in- 
teractions with themselves. That is, we need to make sure that we have the long 
time propagators in a form where they will satisfy the projection operator rule 
GaGa — Ga wherc Go is mixed over orientation. Given that Ga will be a linear 
combination of 1 and au , there are only four possible solutions to this equation 

0, (l±a„)/2, 1. (62) 

These can be found by solving the equation (a -I- PauY — o. + f3au for complex 
constants a and /3. 

To see which of the solutions Eq. is the one obtained when the long 
time propagator Ga — a + jSau interacts with itself we look at the long time 
limit. Noting that ct^ = 1, we find that 

n— ^00 n— ^cxD 2 2 

To make this limit finite and nonzero we renormalize Ga by dividing it by a + /? 
or a — /3, whichever has the larger magnitude. The result will have a limit of 
or 1 only when /3 = 0, a set of measure zero (which does not include our 
propagators). Otherwise the limit will be (1 ± (t„)/2. From this the spin-1/2 
relations given in Eq. ((59)) follow. 

6 Lepton Masses 

In the standard model, mass is an interaction between the left and right handed 
spin-1/2 states. The mass terms in the Lagrangian, for the three generations, 
are: 

^l^limg^l^iPgR + mg^lji^gL). (64) 

where nig is the mass of the gih generation particle. When the ipgL and ipgH are 
rewritten in terms of this paper's long-time propagators, it's natural to expect 
that generational differences in mass will be due to differences in the propaga- 
tors. This provides hope that the three experimentally determined constants, 
rUg can be united into a single constant. The differences between nig will then 
be determined by differences in the corresponding long-time spin propagators. 

The various portions of the long-time propagators differ only in that they 
depend on . Consequently, it may be useful to write rUg in the form: 

mg - E3^iA„«;"9. (65) 
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Since the masses are real, we have that Al = Ai and is real. Putting 
A\ = i?exp(zC) we have: 

mg = A3 + Ay exp(2i5r7r/3) + Ai exp(-2z.g7r/3), , , 

= A3 + 2B cos(257r/3 + C). ^ ' 

In units with c = 1, this is the charged lepton mass equation used by Gerald 
Rosen: 

= 17716 \/eV [1 + \/2cos(2.g7r/3 + 2/9)], (67) 

who notes that it is accurate to O(10~^) and relates it to a Dirac-Goldhaber 
model of the quarks and leptons. [15| In this formula, C = 2/9 and B = 
A-j,^J\j2. The above is similar to the long-time propagators in that it includes 
a square root of 2, but different in that the angle 7r/12 is replaced by 2/9. 
In 1982 3, 17 1, Yoshio Koide discovered a formula for the charged lepton 



masses: 

^2 



2(-yme + -y/TO^ + y/rrir) = 3(me + + m^). (68) 

Since then, measurements of the r mass have converged to Koide's prediction; 
it is still accurate to well within experimental error. This formula follows from 
the 72 in Eq. (|57|. 

The neutrino masses are known by the differences in the squares of their 
masses implied by neutrino oscillation measurements. There are only two such 
measurements, "solar" and "atmospheric", so we need another restriction to 
predict the neutrino masses. Koide's formula would provide the third restriction 
but it is incompatible with the oscillation measurements. On the other hand, 
the oscillation measurements are compatible with Eq. (|66p with A3/B = 
The resulting equation for the neutrino masses 

y/T^ = 0.1000(26)\/cV [1 + \/2 cos(257r/3 + 7r/12 + 2/9]. (69) 

gave predictions for the neutrino masses in 2006: [3] 

m^i = 0.00038, 

m^2 = 0.0089, (70) 
m^3 = 0.0507, 

which give differences in the squares of the masses as: 



VtoL, = \ml,-ml^\ = 7.9 X 10-5 eV^ 



sol 



^"^atm = K2-<3\ - 2.5 X 10-3 eV^ 



(71) 



The mass predictions satisfy Koide's equation but with —^/nTJi. 

A recent measurement of the solar neutrino mass parameter Vm^^j from the 

Sudbury Neutrino Observatory ^ is 7.59(21) x 10"^ eV^ Data from MINOS 
(2O1] give Vm^|_j_^ around 2.43 x 10"'^ eV^, so the above mass predictions are still 
well within the error bars of the oscillation measurements. Adjusting the mass 
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scale 0.1000(26) a little to matcli the new experimental data, a more current 
neutrino mass formula is; 

y/ffh;^ = 0.0990 VeV [1 + \/2cos(2g7r/3 + 7r/12 + 2/9)], (72) 

which gives 7.7 x 10"^ and 2.44 x 10"''' eV^. 

The two lepton mass formulas, Eq. (|571) and Eq. ([7^ . are similar except 
for the mass scale and the angle. The two formulas have the common angle 
2/9, a number which is close to the Cabibbo angle. They differ in that the 
neutrino takes the angle 7r/12 that appears in the long-time projection operators 
Eq. (|49)) . Perhaps the charged lepton mass interaction is a simple interaction 
while the neutrino mass interaction is more complicated. The 7r/12 arises due 
to quantum phase. It may be useful to note that massless spin-1 particles have 
the same quantum phase as spin-1/2 particles, while spin-0 particles have no 
corresponding quantum phase. 

The obvious way to obtain the charged lepton masses in terms of the Ug of 
Eq. (Ull) is: 

^ = 17716 VeV [1 + Uge2^/9-''"/i2 + u^e-^^s+^^/iS] ^ (73) 

the corresponding equation for the neutrinos is: 

= 0.0990 VeV [1 + uge^'/^ + w^e-^^/^], (74) 

The numerical constants in the above two equations, 17716 and 0.0990, have a 
ratio very close to an exact power of three. That is, 17716/0.0990 ~ 311009^ 
speculate that a theory explaining the mass hierarchy between the charged and 
neutral leptons will involve a coupling constant that is a power of 3 (such as the 
1/3 of Eq. (|^5)) ). and that differences in the complexity of the mass interaction 
explains the small neutrino mass. 

7 Discussion 

In ontological theories of quantum mechanics, one attempts to choose which 
mathematical treatment of quantum mechanics is a model of reality in that 
the elements of the model correspond directly to elements of Nature. If we 
choose spin and position as ontological observables, then it is natural for us to 
desire that they have similar relationships between their classical and quantum 
behavior. 

The calculation is robust in that the long-time propagators of Eq. (|49p 
do not depend on the details of Gg. For example, one could consider a theory 
where the paths require all adjacent propagators to be different. This would 
allow paths like XYZX but disallow XXYZ. This change is accomplished by 
putting zeros on the diagonal of Gg. The result would be that G^ and k would 
be different, but the discrete Fourier transform of Gg would still be diagonal 
with one element larger in magnitude than the others and this element would 
dominate G^ leaving G^ unchanged. Spin-1/2 behavior is similarly robust. 
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The calculation uses the spin- up projection operators {X, Y, Z}, and treats 
them equally but does not include the spin-down projection operators {X, Y, Z}. 
Summing over all paths between orthogonal observables such as Z and Z must 
give zero by symmetry considerations. For example, ZXZ+ZXZ = Z{X+X)Z 
= Z{1)Z = 0. Thus one cannot include sums over all paths arbitrarily in a spin 
path integral. More general long-time propagator solutions can be obtained by 
the usual methods of rotation and linear superposition. 

MUBs are a fundamental object of interest in quantum information the- 
ory. They encapsulate the essence of the relationship between complementary 
observables such as position and momentum. To find MUBs at the founda- 
tion of position path integrals, spin-1/2, and the particle generations may be 
more surprising to theorists in elementary particles than to those in quantum 
information. 

In modern elementary particle theory spin arises as a result of examining 



the irreducible representations of the homogeneous Lorentz group [21|, that is, 
spin-1/2 is one of the few possibilities allowed in the intersection of the special 
theory of relativity with quantum mechanics. Physics has found this intersection 
a bountiful place to look for elementary particle models. Despite these successes, 
there has been difficulty combining gravitation with quantum mechanics. 

Our best theories, when extrapolated to very small distances, predict that 
space is anything but flat. Lorentz invariance cannot possibly apply at short dis- 
tances so it cannot logically be used to restrict theories in that regime. Instead, 
the assumptions of this paper must be judged on the basis of how arbitrary they 
are, and whether the resulting calculations are compatible with observations. 

Situations involving extremely non flat spacetime are available, at least theo- 
retically, in that we can look at the dynamics of black holes. Non rotating black 
holes exponentially decay to spherical symmetry. Presumably, their exponential 
approach is accomplished by the radiation of elementary particles. Accordingly, 
L. Motl ^] examined the vibration modes of black holes for various spin cases. 
In addition to the expected results for spin-1/2 and spin-1, he found a spectrum 
of vibrations he called "tripled Pauli statistics" . These were the results of the 
spin-0 and spin-2 vibration modes. Perhaps they have something to do with the 
existence of three generations of elementary fermions. 

8 Acknowledgments 

The author thanks his parents for financial assistance and Marni Sheppeard for 
advice and encouragement, and the anonymous reviewers for improvements in 
the manuscript. 

References 

[1] A. Zee. Quantum Field Theory in a Nutshell. Princeton University Press, 
2003. 



18 



[2] George Svetlichny. Feynman's integral is about mutually unbiased 
bases. Proc. 7th Intl. Conf. Symm. Nonlin. Phys., page 032, jun 2008. 
[quant-ph/0708.3079, 

[3] Jiff Tolar and Goce Chadzitaskos. Feynman's path integral and mutually 



unbiased bases. J. Phys. A, 24:245306, 2009. |quant-ph/0904.0886 



[4] Albert Messiah. Quantum Mechanics. Dover Publications, 1999. 

[5] Byron P. Roe. Particle Physics at the New Millenium. Springer- Ver lag, 
1996. 

[6] Lev Davidovich Landau and Evgeny Mikhailovich Lifshitz. Quantum Me- 
chanics (Non-relativistic Theory). Butterworth-Heinemann, 2002. 

[7] Kaled Dechoum, Humberto de Menezes Franga, and Coraci Pereira Malta. 
Classical aspects of the Pauli-Schrodinger equation. Phys. Lett. A, 248: 
93-102, 1998. labstracti 

[8] Humberto de Menezes Franga. The Stern-Gerlach phenomenon according 
to classical electrodynamics. Found, of Phys., 39:1177-1190. 2009. labstracti 

[9] Brendan Z. Foster and Ted Jacobson. Propagating spinors on a tetrahedral 
spacetime lattice, 2003. ,hep-th/0310166v2, 

[10] Richard P. Feynman and Albert R. Hibbs. Quantum Mechanics and Path 
Integrals. New York: McGraw-Hill, 1965. 

[11] Michael Berry. Quantal phase factors accompanying adiabatic changes. 
Proc. R. Soc. Lond., A392:45-47, 1984. | [abstract] | 

[12] Shivaramakrishnan Pancharatnam. Generalized theory of interference, and 
its applications, part I. coherent pencils. Proc. Ind. Acad. Sci., pages 247- 
262, 1956. 

[13] Rajendra Bhandari. On geometric phase from pure projections. 
J.Mod.Opt, 45(10):2187-2195, 1998. physics/98100 45, 

[14] Tullio Regge. Introduction to complex orbital momenta. // Nuovo Cimento, 
14:951-976, 1959. | [a bstract] | 

[15] Gerald Rosen. Heuristic development of a Dirac-Goldhaber model for lep- 
ton and quark structure. Mod. Phys. Lett. A, 22:283-288, 2007. [[abstract]] 

[16] Yoshio Koide. Fermion-boson two body model of quarks and leptons and 



Cabibbo mixing. Lett. Nuovo Cimento, 34:201, 1982. [abstract 



[17] Yoshio Koide. A fermion-boson composite model of quarks and leptons. 
Phys Lett. B, 120:161-165, 1983. [abstract] . 



[18] Carl Brannen. The lepton masses, 2006. [Link 



19 



[19] Ryan Martin. Results from the neutral current detector phase of the Sud- 
bury Neutrino Observatory, 2009. hep-ex/0905.4907, 

[20] Milind V. Diwan. Recent results from the MINOS experiment, 2009. 
hep-ex/0904.370 6, 

[21] Steven Weinberg. The Quantum Theory of Fields Volume I: Foundations. 
Cambridge University Press, 1995. 

[22] Lubos Motl. An analytical computation of asymptotic Schwarzschild 
quasinormal frequencies. Adv. Theor. Math. Phys., 6:1135-1162, 2003. 
[gr-qc/0212096, 



20 



